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Integrated photonic resonators are widely used to manipulate light propagation in an evanescently-coupled 
waveguide. While the evanescent coupling scheme works well for planar optical systems that are naturally 
waveguide based, many optical applications are free- space based, such as imaging, display, holographies, 
metrology and remote sensing. Here we demonstrate an active dielectric antenna as the interface device that 
allows the large-scale integration capability of silicon photonics to serve the free-space applications. We 
show a novel perturbation-base diffractive coupling scheme that allows a high-Q planer resonator to directly 
interact with and manipulate free-space waves. Using a silicon-based photonic crystal cavity whose 
resonance can be rapidly tuned with a p-i-n junction, a compact spatial light modulator with an extinction 
ratio of 9.5 dB and a modulation speed of 150 MHz is demonstrated. Method to improve the modulation 
speed is discussed. 

I ntegrated photonic resonators are widely used to manipulate light propagation in an evanescently-coupled 
I waveguide, resulting in various on-chip photonic devices including electro-optic modulators 1 ' 2 , filters 3 " 5 , 
I buffers 6 ' 7 and biosensors 8 " 10 . Their compact sizes, low power consumptions and CMOS-compatibility make 
low-cost and large-scale optoelectronic integration possible. The evanescent coupling scheme works well for 
planar optical systems that are naturally waveguide based, such as optical interconnection systems 1112 . But many 
optical applications are free-space based, such as imaging, display, holographies, metrology and remote sensing. 
While the integrated optics field is focusing on building more complex photonic circuits on chip, the planar 
geometry will ultimately limit the capability of these circuits. In the two-dimensional (2D) setting, the massive 
multiplexing capability of optics, which roots from the fact that multiple light beams can propagate in the same 
space without affecting each other, is not taken full advantage of. Extending the dimension from a planar circuit to 
a 3D free-space circuit can significantly expand the capability and scalability of optical systems in data commu- 
nications, processing and storage. 

Here we show a new class of photonic device that can directly manipulate free-space waves. This type of device 
is a bridge for the impressive recent advancement in integrated photonics to be used for free-space optical 
applications and for building compact 3D optical systems. In particular, we demonstrate a silicon-based spatial 
light modulator (SLM) which is the key element in many free-space optical systems 13 . Compared to the existing 
SLMs 14 " 17 , this silicon-based device offers both high speed and the capability of large-scale integration. 

Results 

Diffractive coupling. Conventionally, high-Q optical resonators are designed specifically to prevent light from 
escaping the cavity. Optical coupling to free-space waves is minimized by creating wavevector mismatch between 
the cavity mode and the free-space waves. No matter how well the cavity is designed, however, light can escape to 
free space due to optical scattering caused by side- wall roughness. But the scattering is not strong in any particular 
direction due to the random nature of the roughness. Despite its low intensity, the scattered light can be detected 
from out-of-plane, and the scattering spectrum can be used to determine the resonant wavelength 1819 . Given the 
reciprocity in optics, a free-space wave can couple into the resonator through the reverse process. But the coupling 
efficiency is so low that it can hardly be seen on the transmission or reflection spectrum of the free-space wave. 
Special techniques, such as cross-polarization measurements need to be used to detect the resonance 20 from out- 
of-plane. Obviously, this type of weak coupling cannot be used for SLM. 

We developed a perturbation-based method to convert a photonic resonator to a dielectric antenna that 
strongly couple to free-space waves 21 with coupling efficiency >90%. The design process starts with a conven- 
tional high-Q resonator as the base structure. First, the distribution of an in-plane component of the electric field 
(e.g. E x ) of the resonating mode is obtained from the numerical simulation. Then a perturbation to the base 
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structure is introduced by adding the high-index material where E x is 
positive and removing the high-index material where E x is negative 
(or vice versa). This introduces a perturbation polarization 
Ppert (r) = Ae(r) E(f) that has the same sign over the entire perturbed 
area. In the equation, As is the change in dielectric constant com- 
pared to the based structure. According to the coupled-mode the- 
ory 22,23 , Pp ert at each point can be thought as a source for a secondary 
wave. Since all the secondary waves are in-phase, they interfere con- 
structively and create a strong radiation in the vertical direction. 
Given the reciprocity in optics, the reverse process couples a nor- 
mal-incidence optical beam into the resonator. Here we call this 
coupling method diffractive coupling, which is fundamentally dif- 
ferent from the evanescent coupling method currently used to optic- 
ally access integrated resonators. 

Device design and simulation. The base resonator structure we use 
to construct the SLM is a one- dimensional (ID) photonic crystal 
(PhC) cavity shown in Fig. lb. The ID PhC is formed by periodi- 
cal silicon ribs on a thin silicon slab surrounded by Si0 2 claddings, 
and is designed for the TE mode. The silicon slab is used to inject free 
carriers (electrons and holes) into the cavity from the two sides. A 
high-Q resonating cavity is formed by gradually adjusting the lattice 
constant and thus the bandgap of the PhC 24 " 26 . With a reduced lattice 
constant in the middle of the PhC, as illustrated in Fig. la, the 
bandgap blue-shifts and a confined cavity mode is supported. The 
optical response of the device is simulated by the 3D finite -difference 
time-domain (FDTD) method. The top-view E-field distribution 
(the E x component) of the resonant mode is shown in Fig. lb. One 
can see that £ x has opposite signs in the neighboring ribs. Thus, to 
couple this cavity mode to a normal-incidence optical beam, we 
introduce a perturbation by slightly increasing the widths (in the 
y-direction) of the ribs by Aw in which E x > 0 and slightly de- 
creasing the widths of the ribs by Aw in which E x < 0. The width 
perturbation of Aw is only applied in the cavity area where the lattice 
constant of the PhC is reduced. 

The width perturbation allows the cavity to strongly change the 
transmission and reflection of the normal-incidence waves around 
the resonant wavelength of the resonator. The red line in Fig. lc 
shows the simulated normal-incidence transmission when Aw = 
4 nm, which has a sharp dip with an extinction ratio (ER) of over 



10-dB and a FWHM bandwidth of 0.2 nm (25 GHz). When free 
carriers with the density of 4X10 17 cm" 3 are injected into silicon, 
the reduction in refractive index causes a blue shift of the reso- 
nance 1,31 as shown by the blue line, and a > 10-dB change in the 
optical transmission at the operation wavelength marked by the 
dashed line. 

We simulated the transmission spectra with different perturbation 
strengths Aw, which are shown in Fig. Id. As expected, very little 
resonance feature can be observed without the perturbation (Aw = 
0). As Aw increases, the coupling rate between the cavity mode and 
the normal-incident beam increases. In this process, the cavity Q 
drops since the cavity mode escapes more quickly from the resonator, 
while the ER increases since the coupling rate becomes much higher 
than the intrinsic loss rate. In this set of simulations, the dimensions 
of the devices (see figure caption) are revised to match that of the 
fabricated devices shown below. Given the wafers on hand, the total 
thickness of the silicon layer is dropped to 220 nm, which results in a 
lower off-resonance transmission. In order to make optical align- 
ment easier in the measurements, the fabricated devices are designed 
with a larger cavity area than that in previous simulations. 

Performance determining factors. Compared to the micro -gear 
resonator we demonstrated previously 21 , this PhC cavity provides 
much higher ER and much lower insertion loss due to its sub- 
wavelength periodicity. The off-resonance transmission of the devi- 
ce is higher than 90% (resulting in a low insertion loss) because the 
silicon ribs with sub -wavelength periodicity do not scatter light. To 
the incident light with wavelength far from the resonant wave- 
length, the rib array acts essentially as a uniform effective medium. 
This is in contrast to the micro -gear resonator which acts as a scat- 
tering particle to the off- resonance light due to its wavelength- scale 
dimension. Besides scattering, off-resonance transmission also de- 
pends on layer thicknesses, which will be analyzed in the modeling 
section of this paper. 

The high ER is achieved by eliminating higher-order diffractions 21 
introduced by the perturbation. For the resonating mode of the cavi- 
ty, its far-field radiation pattern is determined by the 2D Fourier 
transform of the perturbation distribution function P pert (x,y), where 
the higher spatial frequency component of Pp ert (x,y) creates radi- 
ation at a larger angle. The spatial frequency component with a 
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Figure 1 | (a) The band diagram of the ID PhC cavity as its spatial period varies along the y-direction. The blue shaded area marks the bandgap. 
(b) Top-view of the PhC cavity where the dashed lines mark the positions of the silicon ribs. The color map shows £ x distribution of the resonant cavity 
mode, (c) The normal-incidence transmission spectra of the resonator without free carriers (red line) and with a free-carrier density 4X 10 17 cm" 3 (blue 
line). The dimensions of the device in (b) and (c) are: the height of silicon slab is 50 nm; the height of silicon ribs is 200 nm; The length of the ribs (in the x 
direction) is 2.5 (am. the lattice constant of the PhC changes gradually from 320 nm outside of the cavity to 280 nm at the center of the cavity in 5 periods. 
Before introducing the perturbation, the rib width is 50% of the local lattice constant. A width perturbation of Aw = 4 nm is introduced inside the cavity 
region. The device is fully surrounded by Si0 2 . (d) The normal-incidence transmission spectra of devices with different perturbations strengths Aw as 
detailed in the legend. The dimensions for this device are: the height of silicon slab is 50 nm; the height of silicon ribs is 170 nm; The length of the ribs is 4 
(am. the lattice constant of the PhC changes from 340 nm outside of the cavity to 290 nm at the center in 7 periods. 
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Figure 2 | Simulated far-field angular distribution of the transmitted light from a micro-gear resonator 21 (a) or the perturbed 1-D PhC cavity (b) when 
they are on resonance. The light-blue ring in (a) shows the l st -order diffraction of the micro-gear resonator, which is absent from the radiation of the PhC 
cavity. 



period shorter than the wavelength do not create far-field radiation, 
therefore the radiation pattern is determined by the slow-varying 
profile of Pp ert (x,y). 

In the micro -gear structure, the perturbation locates at the inner 
and outer edges of the rings with radii comparable to the wavelength. 
Its far- field radiation pattern obtained from 3D FDTD simulation is 
shown in Fig. 2(a). Besides the O th -order diffraction that is used to 
couple the normal-incidence light, there is also a l st -order diffraction 
that acts as an addition source of loss for the cavity which limits the 
extinction ratio of the resonance. In contrast, P pert (x,y)m this PhC 
cavity distributes at the edges of the ribs with a deeply sub-wave- 
length periodicity. Therefore no higher-order diffraction is expected, 
which is confirmed by the far-field radiation pattern shown in 
Fig. 2(b). The optical Q of the PhC cavity is optimized when the 
cavity mode has a Gaussian profile 27 . Under that condition, P pert 
introduced by a constant Aw also has a Gaussian profile and the 
far-field radiation pattern from the cavity matches well with that of 
a normal-incident Gaussian beam. This would maximize the coup- 
ling efficiency to the resonating mode and the extinction ratio on the 
transmission spectrum. If coupling to other types of beams is desired, 
a different distribution of Aw can be used. 

Measured spectra from fabricated devices. SLMs with different 
perturbation strengths Aw are fabricated on a silicon- on-insulator 
(SOI) substrate by deep-UV lithography in a CMOS photonics 
foundry (See Methods for details). As shown in Fig. 3a, a p-i-n 
junction is built across the PhC cavity to control the carrier density 
in the cavity. Fig. 3b shows an SEM picture of the device taken after 
the Si0 2 top cladding layer is removed by HF wet etching. The design 
parameters for the fabricated devices are: the height of silicon slab is 
50 nm; the height of silicon ribs is 170 nm; The length of the ribs (in 
the x direction) is 4 |im. The lattice constant of the PhC changes 
gradually from 340 nm outside the cavity region to 290 nm at the 
center of the cavity in 7 periods. Without perturbation, the rib width 
is 50% of the local lattice constant. 

The measured normal-incidence transmission and reflection 
spectra of devices with different width perturbation Aw = 0 nm, 
2 nm and 4 nm are shown in Fig. 4a and 4b respectively. As we 
expected, no resonance feature can be identified if no perturbation 
is introduced. On the other hand, even a perturbation as small as 
2 nm can introduce a strong resonant feature that is clearly observed. 
The transmission drops over 90% at the resonant wavelength, show- 
ing that over 90% of optical power has been coupled into the res- 
onator, which then escape through scattering or reflection. The 
device with a larger perturbation shows a stronger coupling, a 
broader bandwidth and a higher ER. The resonant spectra have an 
asymmetric shape due to the interference from the partial reflections 



at the multiple dielectric interfaces through the SOI substrate, which 
is not taken into consideration in the FDTD simulation. 

Fitting to theoretical model. The measured spectra agree well with 
those calculated from an analytical model based on the coupled- 
mode theory and the transfer matrix method. The device built on 
the SOI substrate can be modeled as a multilayer dielectric structure 28 
shown in Fig. 5a. The /-th dielectric interface is described by a 
(«/_ i + m)/2ni (m - m-i)/ 2ni 
_ (m - ni- i)/2n/ + »/)/ 2nj 

th layer is described by a transfer matrix (also called propagation 

-ik 0 nid 



transfer matrix T/ = 



matrix) P/ = 



and the /- 



^ _ mmiCl , where n t is the refractive index of 

the /-th layer and d\ is the thickness of the /-th layer. The relationship 
between the amplitudes of the optical waves traveling downward and 



upward in the (M)-th layer and the /-th layer is thus 



= P V T V 



E + 
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where the and the ' signs mark the 



downward and upward travelling waves respectively. To model the 
non-resonant behavior of the silicon rib array, the rib layer is 
modeled as a homogeneous layer with an effective index n eff that 
depends on the occupation ratio of the silicon ribs and the 
polarization of light. 

The resonance of the PhC cavity is modeled as a virtual interface in 
the middle of the rib array layer, and is described by a transfer matrix 
T c , which can be obtained from the coupled-mode theory. Due to the 
perturbation induced diffraction, the PhC cavity mode couples to 
both the upward and downward traveling optical beams. It can be 



p + -doped 




Figure 3 | (a) The diagram of the SLM design where a p-i-n junction is 
built across the ID PhC cavity, (b) A top-view SEM picture of the 
fabricated device after the Si0 2 upper- cladding is removed by wet etching 
in HF. 
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Figure 4 | Measured normal-incidence transmission and reflection spectra of the SLMs with different perturbation strengths Aw. (a) The transmission 
spectra, (b) The reflection spectrum. In both panels, the symbols are the measured data with Aw = 0 nm (black dots), Aw = 2 nm (red circles) and Aw = 
4 nm (blue dimonds). The solid lines are the fitted spectra using the theoretical model with Aw = 2 nm (red line) and Aw = 4 nm (blue line). The spectra 
are normalized to the transmission and reflection when the optical beam passes through a featureless area of the chip. The featureless area has a 
transmission of —50% compared to air (i.e. when the chip is completely removed from the setup), which is mainly due to the reflection from the silicon 
substrate and optical scattering at the unpolished backside of the chip. 



express as a single-mode optical resonator coupled with two ports as 
shown in the inset figure of Fig. 5a. The dynamic equations for the 
amplitude of the resonance mode can be written as 29,30 

|2n 



da 
di 



i(oj-o)o)- 



1 

2t 0 



I*i| 2 + 



K 2 \ 



CI ~~\~ KjS- 



2*2 



-K,a 



(i) 

(2) 
(3) 



where w is the frequency of the input wave and w 0 is the resonant 
frequency of the PhC cavity. In these equations, s describes the 
amplitude of the input or output waves, which is normalized such 
that \s\ 2 is the power of the wave, a describes the amplitude of the 
resonant mode, which is normalized such that \ a\ 2 is the energy in the 
cavity. t 0 is the intrinsic photonic lifetime of the cavity, which deter- 
mines the intrinsic quality factor as t 0 = Qiko/lnc. K^and k 2 are 
the field coupling coefficients associate with the downward (s^ ,s/ ) 
and upward (sf ,s^~ ) propagating waves. The coupling quality factor 
Q c is determined by the coupling coefficients as 



|Kl| 2 + 



k 2 \ =2to/Q c A 0 



(4) 



Under the steady state (da/dt = 0), the amplitude a of the resonant 
mode can be expressed as 

_ (kisf +k 2 s~) 



- i(co - wo) + ( | h | 2 + | k 2 1 2 ) /2 + 1 /2t 0 



(5) 



By plug Eq. (5) into Eq. (2) and (3), one can obtain the transfer matrix 



of the resonator T c as 



= T r 



As the propagation property of each layer and their interfaces are 
described by a (2X2) transfer matrix, the behavior of the entire 
multilayer structure is determined by the product of these matrices, 
from which the transmission and reflection spectra of the device are 
calculated. The ER of the resonance typically increases with the ratio 
Qi/ Q c . The basic parameters of the resonator are obtained by fitting 
the calculated spectra to the measured ones. For the fabricated device 
with Aw = 2 nm, the coupling strength \k\ |and \ k\ \ are found to be 
1.02X10 11 s" 1 and 1.17X10 11 s" 1 respectively, which gives a coupling 
Q (Q c ) of —5,500. The intrinsic lifetimeTo is 16.4 ps, corresponding to 
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Figure 5 | (a) A diagram of the multilayer dielectric model of the SLM. The multi-layer structure is described by the transfer matrics P/ and T b where 
/= 1,2. .6. The resonant cavity is described by the transfer matrix T c . (b) Off-resonance transmission versus the total thickness of the silicon layer and the 
thickness of the slab layer. The diamond and triangle symbols marks the thicknesses used in the FDTD simulations in Fig.lc and Fig. Id, respectively. 
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Figure 6 | Measured normal- incidence transmission (a) and reflection (b) spectra of the SLM with a perturbation of Aw = 2 nm under different bias 
voltages V DC . Blue line: V DC = 0-7 V; green line: V DC = 1.0 V; red line: V DC = 1.1 V. The solid lines are measured when the E-field of the normal- 
incidence light is parallel to the silicon slabs (x-polarized). The dashed black line is measured when the E-field is perpendicular to the silicon slabs (y- 
polarized). The spectra are normalized to the transmission and reflection when the optical beam passes through a featureless area of the chip. 



an intrinsic quality factor Q f of —20,000. The coupled photon life 
time in the cavity is 3.6 ps, corresponding to a total Q of Q = 
1/(1/Qj + 1/Q c )~4,300. When the width perturbation increases 
to be 4 nm, Q z and Q c decrease to 7,000 and 2,300 respectively. 

Based on this model, the off-resonance transmission for different 
thickness of the total silicon layer and slab layer is calculated, which is 
shown in Fig. 5b. The two marked points represent the two cases that 
are used in the simulations for Fig. lc and Fig. Id, respectively. One 
can see that the off-resonance transmission can be higher than 99% 
which optimal layer thicknesses. 

High-speed modulation. When a forward bias voltage is applied on 
the p-i-n junction, free carriers are injected into the cavity which 
changes the refractive index of silicon and blue -shift the reso- 
nance 1 ' 31 . For a device with Aw = 2 nm, the measured normal- 
incidence transmission and reflection spectra at different DC bias 
voltage are shown in Fig. 6a and Fig. 6b, respectively. With a low 
carrier density (V DC — 0.7 V), the transmission spectrum is the same 
as that at V DC = 0 V, which shows a resonant dip with an ER of 
9.5 dB. When free carriers are injected into the p-i-n junction with a 
forward bias voltage of 1.0 V, the resonance blue-shifts. Fitting of the 
spectra gives a decreased intrinsic Q of 1 1,000 due to the free-carrier 
absorption effect, while the coupling Q does not change as one 
expects. The extinction ratio of the resonance drops due to the 
additional loss in the cavity from the free-carrier absorption effect. 
The resonances can be seen only when the input light is x-polarized 
(E-field parallel to the ribs), which interacts with the TE mode of the 
cavity. The y-polarized light does not couple to the resonant mode, 



therefore its transmission and reflection are nearly constant, as the 
dashed lines in Fig. 6a and 6b show. The transmission and reflection 
of the jy-polarized light are different from the off- resonance 
transmission and reflection of the x-polarized light, because the 
effective index of the sub -wavelength grating as an effective 
medium is different for the two polarizations. 

Due to the asymmetrical line shape of the resonant spectra, the 
transition from low to high transmission happens in a wavelength 
range much narrower than the total Q of the resonator suggests. Here 
we define an effective Q to describe the sharpness of the transition 
edge, which determines the modulation efficiency of the SLM. In the 
measured device with Aw = 2 nm, Q e ff = A 0 /A/L~ 7,500, where k Q is 
the resonant wavelength and Ak is the wavelength different between 
the minimal and maximal transmission points, as shown in Fig. 6a. 

When the laser wavelength is fixed at 1547.4 nm and the bias 
voltage on the p-i-n junction is scanned, the transfer function is 
plotted in Fig. 7a. Transmission changes significantly when the bias 
voltage changes only 0.1 V. To show the high-speed modulation, a 
150-MHz pseudo-random data sequence is applied on the p-i-n 
junction and the waveform of the transmitted optical signal is shown 
in Fig. 7b. 

Discussion 

The speed of the SLM is limited by the width of the p-i-n junction in 
this device where the intrinsic region is ~5-jim wide. The relatively 
wide intrinsic region limit the strength of the electric field inside 
the junction for carrier exaction, and electrons have to travel a longer 
distance to go out of the wider junction. Since the optical field 
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Bias Voltage V DC (V) Time(ns) 

Figure 7 | Electro-optic performance of the SLM. (a) Normalized optical transmission of the SLM at the wavelength of 1547.4 nm versus the bias voltage 
applied on the p-i-n junction, (b) The waveform of the modulated optical output at the bit-rate of 150 Mbit/s. 
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Table 1 


| Coupling Q (Q c ) 


and intrinsic Q (Q,) of the simulated and 


measured devices with different perturbation strengths Aw 






Simulated 




Measured 


Aw 
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2 nm 
30,860 
1 70,000 


4 nm 8 nm 
7,346 1,837 
1 00,000 60,000 


12 nm 
994 
55,000 


2 nm 4 nm 8 nm 
5,500 2,300 810 
20,000 7,000 1 ,000 



intensity is minimal in the areas between the silicon ribs, as shown in 
Fig. lb, much narrower p-i-n junctions can be built without signifi- 
cantly increasing the cavity loss by extending the p- and n- doped 
regions in to these areas in an interleaving fashion. Multi-GHz 
modulation 32 can be achieved with such a design. 

The measured and simulated spectra of devices with different 
perturbation strengths are fitted with the theoretical model, from 
which the intrinsic Q (Q f ) and coupling Q (Q c ) are extracted and 
listed in Table 1. From the simulation results, one can see that the Q c 
roughly scales inversely with the square of perturbation strength (i.e. 
Q c ~ Aw -2 ). This is consistent with the prediction of Eq. (4) where 
the field coupling efficiencies K^and k 2 are expected to increase lin- 
early with Aw. The intrinsic Q of the fabricated device is lower than 
that of the simulated one. This is mainly due to the scattering loss 
from the side-wall roughness and other fabrication imperfection. 

Besides that, there are two phenomena that are not expected from 
the basic theory. One is that Q, drops as the perturbation increases, 
which can be seen from both simulations and the experiments. We 
believe this drop is just an artifact from the fitting process. In the 
fitting process, the ratio Q//Q c is adjusted to have the ER of the 
calculated spectra match that of the simulated or measured ones. 
In theory, ER increase with the ratio Q//Q c . However, ER in the 
simulated and measured spectra also depends on how well the far- 
field radiation pattern of the cavity matches the spatial profile of the 
input optical beam. The drop of ER due to an imperfect match caps 
Qi/Q c in the fitting process, so that Q z appears to decrease with Q c . 
The fabricated device with Aw = 8 nm shows a more significant drop 
in Qi. That is because the width of the narrowest trenches in this 
structure falls just below the resolution of the optical lithography 
process, so that some trenches are not completely cleared, which 
introduces extra optical loss in the cavity. The other unexpected 
phenomenon is that the diffractive coupling in the measured device 
is much stronger (represented by a lower Q c ) than that obtained in 
the simulation with the same perturbation strength. It is yet to know 
whether this is due to a systematical under-estimation of the coupling 
strength in the simulation, or due to that the fabricated device has a 
larger perturbation than what is designed. 

Methods 

SLMs are fabricated in a CMOS photonics foundry at the Institute of Microelectronics 
of Singapore. The fabrication starts on an SOI wafer with a 220-nm-thick silicon layer 
and a 3 -urn-thick buried oxide layer. To construct the ID PhC cavities, silicon ribs 
with the height of 170 nm are patterned on a silicon slab with the thickness of 50 nm 
using 248-nm deep-UV lithography and inductively- coupled plasma etching. 
Following the etching, the p-i-n junctions are formed by patterned ion implantations 
with a dosage of 5X 10 14 cm" 2 for both the p + and n + doping regions. A 2.1-|im-thick 
Si0 2 layer is then deposited onto the wafer using plasma-enhanced chemical vapor 
deposition (PECVD). Finally, vias are opened on the ion-implanted areas and a 1.5- 
um-thick aluminum layer is sputtered and etched to form the electric connections. 
The serial resistance of the diode is measured to be 105 Q. After the fabrication 
process, the contact pads connecting to the p-i-n junction are wire-bonded to a SMA 
connector with a 50-ohms terminal resistor for impedance match. 

The 3D FDTD simulations are done with commercial software Lumermical FDTD. 
A non-uniform grid is used which has a spatial resolution —30 nm around the 
resonator. Even though perturbation we introduced is much smaller than the grid 
size, the software is capable of incorporate that in the simulation. When a dielectric 
interface (Si/Si02) lies between two grid points, the program modifies the dielectric 
constant at the neighboring grid points according to the position of interface. This 
way, the small shift of the dielectric interface due to the width perturbation is taken 
into account in the simulation. 

The transmission and reflection spectra are measured with a tunable laser whose 
fiber-coupled output is collimated to form a free-space optical beam. The polarization 



of the light is adjusted by a polarization controller. Then the optical beam is focused 
onto the device from top with an aspheric lens. The transmission through the back of 
the chip is collimated by a lens and measured with a photodetector. The reflection 
spectra are taken by inserting a beam splitter in the path of the incident light beam. 
Electro-optic modulation is achieved by applying DC and RF voltages on the p-i-n 
junctions through the SMA connector. 
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